Refractive effects in the scattering of loosely bound nuclei 
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A study of the interaction of loosely bound nuclei ^''^Li at 9 and 19 A MeV with light targets 
has been undertaken. With the determination of unambiguous optical potentials in mind, elastic 
data for four projectile-target combinations and one neutron transfer reaction ^■^C(^Li,*Li)^^C have 
been measured on a large angular range. The kinematical regime encompasses a region where 
the mean field (optical potential) has a marked variation with mass and energy, but turns out 
to be sufficiently surface transparent to allow strong refractive effects to be manifested in elastic 
scattering data at intermediate angles. The identified exotic feature, a "plateau" in the angular 
distributions at intermediate angles, is fully confirmed in four reaction channels and interpreted as 
a pre-rainbow oscillation resulting from the interference of the barrier and internal barrier far-side 
scattering subamplitudes. 

PACS numbers: PACS number(s): 25.70.Bc, 24.10.Ht, 25. 70. Hi, 27.20.+n 
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I. INTRODUCTION 

The study of nucleus-nucleus elastic scattering has a 
long history and remains of interest due to both suc- 
cesses and failures that mark it (see for example, Refs. 

2] and references therein). It is an important subject 
per se, and is also important as a tool for the descrip- 
tion of a series of phenomena that involve the distorted 
waves given by optical model potentials (OMP). We are 
searching here for reliable ways to predict optical model 
potentials for reactions with radioactive nuclear beams 
(RNB). In particular our interest focuses on finding re- 
liable descriptions for transfer reactions involving rela- 
tively light, loosely bound nuclei, which are used in in- 
direct methods in nuclear astrophysics. A range of RNB 
studies were made at energies around 10 MeV/nucleon, 
where the reactions are peripheral, with the intent to ob- 
tain information about the surface of the nuclei involved. 
These reactions use DWBA techniques to extract nuclear 
structure information. However, the well known exis- 
tence of many ambiguities in the OMPs extracted from 
elastic scattering can raise questions about the accuracy 
of these determinations. Therefore, we are searching for 
ways to reduce these ambiguities and to predict OMP for 
reactions with RNBs. Experimental studies using RNBs 
have, heretofore, not been suitable for detailed elastic 
scattering analyses. The closest we can get using sta- 
ble beams is by studying the elastic scattering of loosely 
bound nuclei. We chose here to study the elastic scatter- 
ing of ^'"^Li projectiles, because they are fragile (loosely 
bound), with a pronounced cluster structure and with 
low Z and can, therefore, exhibit a range of phenomena, 
involving absorption, diffraction and refraction, mostly 
of nuclear nature. 

Earlier we have carried out a study of elastic scattering 
around lOA MeV for a range of projectile-target combi- 



nations involving p-shell nuclei . We found a relatively 
simple method to predict OMP for loosely bound nuclei, 
based on the renormalization of the independent real and 
imaginary terms obtained from a double folding proce- 
dure using the JLM nucleon-nucleon (NN) effective inter- 
action. The procedure successfully described the data for 
all the projectile-target combinations and the energies in 
the study for most of the angular ranges measured. In 
one single case CLi at 99 and 130 MeV on ^^C) the fold- 
ing potentials failed to describe well the large angle data. 
Later, the results were used to describe elastic scattering 
angular distributions measured in a series of experiments 
with RNBs at or around 10 MeV/nucleon: ''Be on ^°B 
and melamine targets 0, "C 0, [| and ^^F 
on ^^C and ^^N targets. We return to that study here 
with new data extending the angular ranges for the ^Li 
scattering and adding data for ^Li scattering and with a 
refined analysis. 

Recent work "3, has established that elastic scat- 
tering of light tightly bound heavy ion systems such as 
^^0-f ^■^C and ^^O-f-^^O show sufficient transparency for 
the cross section to be dominated by the far-side scatter- 
ing. Intermediate angle structures appearing in the elas- 
tic scattering distributions at angles beyond the Fraun- 
hofer diffractive region have been identified as Airy min- 
ima of a nuclear rainbow, i.e. a destructive interference 
between two far-side trajectories which sample the inte- 
rior of the potential. A number of high order Airy min- 
ima have been identified by observing that such struc- 
tures are largely insensitive to an artificial reduction of 
the absorption in the optical potential, and therefore they 
appear as a manifestation of the refractive power of the 
nuclear potential. While at high energy (T^ this picture 
was well substantiated by a semiclassical nonuniform de- 
composition of the scattering function |ll| , at lower ener- 
gies the situation is more difficult to understand. It has 
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been shown by Anni that such structures could be 
explained by the interference of two amplitudes appear- 
ing in different terms of a multireflection uniform series 
expansion of the scattering amplitude and therefore the 
interpretation using rainbow terminology is not appro- 
priate. 

For loosely bound nuclei the situation is even more un- 
certain. When a nucleon or a group of nucleons has small 
separation energy, the wave function penetrates well be- 
yond the potential range. The corresponding components 
in the optical potential are expected to be more diffuse as 
compared to normal nuclei, leading to a competition be- 
tween the increased refractive power of the real potential 
and the increased absorption at the nuclear surface. The 
small separation energy implies also that the dynamic 
polarization potential (DPP) ^l^J arising from the cou- 
pling to breakup states may be strong and with a com- 
plicated energy and radial dependence. It follows that 
for loosely bound nuclei the DPP cannot be treated as a 
small perturbation and the usual phenomenological pro- 
cedure in renormalizing the folding potential form factor 
may be questioned. It has been estimated that the DPP 
is strongly repulsive at the nuclear surface in the case 
of ^Li [iJI and this prompted Mahaux, Ngo and Satch- 
ler ^5 to conjecture that for loosely bound nuclei the 
barrier anomaly may be absent due to the cancellation 
between the repulsive (DPP) and attractive (dispersive) 
components of the optical potential. 

In the specific case of ^'^Li scattering on light targets a 
large body of data have been accumulated in the range 5- 
50 MeV/nucleon. At high energy, Nadasen and his group 
|16, 17] have been able to derive a unique optical poten- 
tial which was essential to assess the quality of the folding 
model. At lower energies, ambiguities found in the anal- 
ysis of data prevented any definite conclusion about the 
strength and energy dependence of the optical potential. 
A study by Trcka et al.^^ on ^Li-f^^C elastic scatter- 
ing at 50 MeV, found an exotic feature ("plateau") in 
the angular distribution of the elastic scattering at inter- 
mediate angles which resembles similar structures found 
in more bound systems. They interpreted the structure 
as a diffractive effect arising from an angular momen- 
tum dependent absorption. There are experimental hints 
that such structures also appear in neighboring systems, 
^Li-|-^^0 and ^Li-f ^Be, as a possible manifestation of the 
average properties of the interaction potential. 

In this paper we present a precision measurement of 
elastic scattering of ^^^Li on 12,13q ^^^^^ 9gg targets at 
9 and 19 MeV/nucleon. The lower energy was chosen 
in view of our systematic studies of nuclear reactions for 
astrophysics. The higher energy is close to the satura- 
tion energy for these projectiles, i.e. the energy where 
almost all reaction channels are open. The "plateau" 
feature is confirmed in four projectile-target combina- 
tions at 9 MeV/nucleon. The high selectivity induced by 
this structure allowed the derivation of an almost unique 
Woods-Saxon optical potential. A folding model analysis 
using the complex, density and energy dependent NN in- 



teraction of Jeukenne, Lejeune and Mahaux (JLM) [l^ . 
where corrections due to the strong DPP have been in- 
cluded, confirmed that our elastic distributions could be 
described using deep and extremely transparent poten- 
tials. The remaining ambiguities have been eliminated 
using an accurate dispersion relation analysis. The in- 
termediate angle structures have been discussed using 
the semiclassical uniform approximation for the scatter- 
ing function of Brink and Takigawa . We explain the 
intermediate angle structure as a coherent interference 
effect of two subamplitudes corresponding to trajectories 
refiected at the barrier and interfering with trajectories 
which sample the nuclear interior. Thus, this refractive 
effect appears as a signature of a highly transparent in- 
teraction potential. 

The paper is structured in the following way: after this 
introduction, the experimental methods are discussed in 
Sect, m the analysis of the elastic scattering data using 
phenomenological and microscopic optical model poten- 
tials is discussed in Sect. IIIII and the implications of this 
analysis for the transfer reaction (^Li,*Li) is discussed 
in Sect. IIVI In Sect. 0the dispersion relation is used 
to put additional constrains on the potentials extracted, 
followed by a discussion of the decomposition of the far- 
side scattering amplitude into barrier and internal barrier 
components responsible for the "plateau" structure at in- 
termediate angles (Sect. IVI|I . and the conclusions (Sect. 

Ivnl) . 



II. THE EXPERIMENTS 

The experiments were performed using ®Li and ^Li 
beams of 9 and 19 MeV/nucleon from the Texas A&M 
University K500 superconducting cyclotron and the Mul- 
tipole Dipole Multipole (MDM) magnetic spectrometer 
|2ll | . A list of the measurements is given in Table |l| The 
measurements with ''Li were done to extend the angular 
range that was effectively covered in earlier work. The 
experimental setup and the data reduction procedures 
were similar to those used in Ref. 0. The beams were 
prepared using the beam analysis system f"2.2j , which al- 
lows for the control of the energy spread {AE/E up to 
1/2500) and angular spread (0.1°) of the beam. Self- 
supported ^Be (200 fig/cnP thick), (260 fig/cm^) 
and (390 ^g/cm^) targets were placed perpendic- 
ular to the beam in the target chamber of the MDM. 
The magnetic field of the MDM spectrometer was set to 
transport the fully stripped Li ions to the focal plane 
where they were observed in the modified Oxford detec- 
tor |2^. In the detector, the position of the particles 
along the dispersive direction was measured with resis- 
tive wires at four different depths, separated by about 
16 cm each. For particle identification we used the spe- 
cific energy loss measured in the ionization chamber and 
the residual energy measured in a NE102A plastic scintil- 
lator located behind the output window of the detector. 
The input and output windows of the detector were made 
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of 1.8 and 7.2 mg/cm^ thick Kapton foils, respectively. 
The ionization chamber was filled with pure isobutane 
at 40 torr. The entire horizontal acceptance of the spec- 
trometer, /S.9 = ±2°, and a restricted vertical opening, 
A0 — ±0.5°, were used in the measurements at forward 
angles, whereas at the largest angles the vertical open- 
ing of the acceptance window was raised to = ±1.0°. 
Raytracing was used to reconstruct the scattering angle. 
For this purpose, position calibration of the detector was 
performed by using the scattering on a thin Au target 
(212 /ig/cm^) and an angle mask consisting of five open- 
ings of 59^0.1°, located at -1.6°, -0.8°, 0°, +0.8° and 
+1.6° relative to the central angle of the spectrometer. 
In addition to RAYTRACE Q calculations, angle cal- 
ibration data were obtained at several angles by using 
the angle mask. Typically the spectrometer was moved 
by 2° or 3° at a time, allowing for an angle overlap that 
provided a self-consistency check of the data. Normal- 
ization of the data was done using current integration 
in a Faraday cup. Focal plane reconstruction was done 
at each angle using the position measured with the sig- 
nals in the wire nearest to the focal plane and using the 
detector angle obtained from the position measured at 
two of the four wires (typically the first and last). The 
angular range, — 4°, covered by the acceptance slit 
was divided into 8 bins, resulting in 8 points in the an- 
gular distribution being measured simultaneously, with 
each integrating over 69iab — 0.5°. 

The measurements with the angle mask showed that 
the resolution in the scattering angle (laboratory) was 
Mres = 0.18° - 0.25° full- width at half maximum 
(FWHM) . This includes a contribution from the angular 
spread of the beam of about 0.1°. The best energy reso- 
lution obtained at forward angles was 150 keV FWHM. 
It degraded as we advanced to larger angles due to the 
large kinematic factor, k — coupled with the finite 
angular spread in the beam. However, it never degraded 
our ability to isolate the elastic peak, even in the case 
of the ^Li experiments where the first excited state of 
the projectile is only 477 keV away. The active length 
of the focal plane allowed us to cover a total excitation 
energy of about 7 MeV, centered around the elastic peak. 
Thus we were able to measure inelastic scattering to the 
lowest excited states of the projectile-target systems at 
the same time. These inelastic scattering data were used 
as additional information to check the experimental pro- 
cedures. In one of the runs we have also measured the 
neutron transfer reaction '^^C{^lj\,^LiY'^C at E('^Li)=63 
MeV, which was discussed elsewhere in detail and is 
used here to check the sensitivity of observables in other 
channels to the OMP extracted from the elastic scatter- 
ing data. 

To obtain accurate absolute values for the cross sec- 
tions, target thickness and charge collection factors were 
determined by a two-target method as described in Ref. 
[2^ . We also determined the target thickness by measur- 
ing the energy loss of alpha particles from a ^^*Th source 
and the accuracy in normalization is 9%. Combining the 



results of these independent determinations, we conclude 
that we have an overall normalization accuracy of 7% for 
the absolute values of the cross sections. 



III. OPTICAL-MODEL ANALYSIS 

The measured elastic data at 9 MeV/nucleon, shown 
in Fig. ^as the ratio to the Rutherford cross section, ex- 
tend to a larger angular range than previously measured 
Q. These data show complex forms with characteristic 
rapid oscillations at small angles followed by a marked 
change in shape at intermediate angles: a plateau is de- 
veloped aX 9 = 50° — 70° followed by a deep minimum 
ai 9 K, 80°. Assuming pure Fraunhofer scattering at for- 
ward angles, we extract a grazing angular momentum 
Ig «15 from the angular spacing /S.9 — i^jilg + 1/2). 
The striking fact is that the same pattern emerges for 
all four projectile-target combinations, including that for 
the ^Be target where a much stronger absorption is ex- 
pected. (We remind the reader that ^Be is a perfect black 
disc target since it has very low thresholds for breakup 
into the neutron and alpha channels S'„=1.66 MeV and 
5q.=2.47 MeV, and there are no bound excited states. 
These values should be compared with 5q,=1.47 MeV in 
^Li and 5*0=2.47 MeV in ^Li.) 

Similarities seen in the differential cross sections shown 
in Fig. n indicate general wave-mechanical character- 
istics of the scattering process and average systematic 
properties of the nuclear interaction. Specific structure 
effects can be isolated only as small deviations from the 
normal behavior. Therefore the data are analyzed using 
optical potentials with conventional Woods-Saxon (WS) 
form factors for the nuclear term, supplemented with a 
Coulomb potential generated by a uniform charge dis- 
tribution with a reduced radius fixed to r^—^ fm. No 
preference has been found for volume or surface localized 
absorption and throughout the paper only volume ab- 
sorption is considered. In the absence of any spin depen- 
dent observables, spin-orbit or tensor interactions have 
been ignored. Ground state reorientation couplings also 
have been neglected. The potential is defined by six pa- 
rameters specifying the depth and geometry of the real 
and imaginary terms 



U(r) = -(Vfv(r) + iWfw(r)) 



where 



fx(r) 



1 + exp 



'r-r,(A;/-VA^/3) 



(1) 



(2) 



and x=V,W stands for the real and imaginary parts of 
the potentials, respectively. The number of data points 
per angular distribution exceeds N=100 points and there- 
fore the usual goodness of fit criteria (x^) normalized to 
N has been used. A source of bias was the finite angular 
acceptance of the detectors (the 0.5° bins, in the present 
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FIG. 1: (Color online) Woods-Saxon optical model analysis (full lines) of elastic scattering data (open points) at 9 MeV/nucleon 
(Table Far-side/near-side cross sections are also shown by dashed and dotted lines, respectively. The depth of the real 
potential is shown to identify the particular WS potential parameters used in the calculations. 



case). The averaging associated with this finite angular 
resolutions has most effect on the depth of sharp minima. 
A few exploratory calculations showed that allowing the 
normalization to vary did not result in any qualitative 
changes and did not indicate that any renormalization 
by more than a few percent would be preferred. Optical 
parameter sets collected from the literature were used as 
starting values for the search procedure. In particular the 
potential OMl of Trcka et al. 18] has been extensively 
tested. Guided by these potentials and by our earlier 
analysis Q a number of some 10^ potentials with real 
volume integrals in the range Jv = 200 — 600 MeV fm^ 
have been generated for each reaction channel, thus ex- 
ploring the functional Woods-Saxon space in full detail. 
Local minima were identified and a complete search on 
all six parameters determined the best fit potentials. The 



plateau feature at the intermediate angles and the sharp 
decrease in the cross section near 9 = 80° could be fit- 
ted only with deep potentials with real volume integrals 
(per nucleon) exceeding a critical value Jvcrit ~ 300 MeV 
fm'^. There is a consistent preference for potentials with 
relatively weak imaginary parts with values of W around 
15 MeV except for ^Li scattering where somewhat larger 
values are needed to fit the data. We systematically find 
rv < rw and large diffuseness parameters ay — cLw — 0.8 
fm in agreement with theoretical expectations for loosely 
bound nuclei [23, HI]. A grid search procedure on the 
real depth of the potential allowed to identify discrete 
ambiguities. Parameters for the first two discrete fami- 
lies are given in Table UTI These are identified by a jump 
of A Jy « 100 MeV fm"^ from one family to another and 
almost constant imaginary volume integral. As a con- 
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sequence, the total reaction cross section seems to be a 
well determined observable. Gridding on other WS pa- 
rameters revealed a continuous ambiguity of the form 
JyRv ~ const, where Ry is the rms radius of the po- 
tential. The larger the volume integral, the smaller the 
radius that is required to fit the data. This is a clear 
manifestation of a complicated radial dependence of the 
dynamic polarization potential (DPP) which may lead to 
radii much smaller than the minimal value implied by the 
folding model (e.g., Rj^ = Rf + R2, for a zero range NN 
effective interaction). However, for each discrete family 
rather precise values of the rms radii were required to fit 
both forward and intermediate angle cross sections. 

Sometimes more subjective criteria may be used to 
choose between various ambiguous potentials based upon 
general theoretical expectations. For example one may 
require consistency with the results of analyses of other 
data for the same system at nearby energies with the 
expectation that the potential should not change rapidly 
with mass and energy. Individual elastic data sets possess 
individual idiosyncrasies which facilitate the inference of 
a single local potential. We note that, seemingly, there is 
a compatibility between all data sets: an optimum poten- 
tial found for one data set gives already a good fit to the 
other. In fact, potentials given as first entry in Table UTI 
were obtained by iterating several times this procedure in 
an attempt to find a single potential which would simul- 
taneously fit all data at 9 MeV/nuclcon. A compromise 
could be obtained with transparent deep potentials close 
to Vb ~ 225 MeV having a strongly refractive core at 
small radii, surrounded by a weakly absorptive halo. In 
fact, examining the ratio w{r) = W{r)/V{r) as a 
function of the radial distance, we found that our poten- 
tial may by qualified as having internal (r ^ — 4 fm) 
and surface (r > 8 fm) transparency {w ~ 0.1) but with 
a pronounced maximum {w ~ 0.8) near the empirical 
strong absorption radius « 6 fm) in agreement with 
the systematics found in other more bound systems |29j |. 
The surface localized absorption suggests that the reac- 
tion mechanism is dominated by direct reactions. The 
relatively large radius of the absorption required by the 
data is an indication that fusion already sets in the re- 
gion of the barrier and that fusion is a large component 
of the total reaction cross section. 

A variety of notch tests have been performed to deter- 
mine the radial sensitivity of the potential. One test was 
done using a Gaussian spike superimposed on the real 
potential at a given radius. The resulting influence on 
the of the fits is displayed in Fig. |2J It shows that 
there is a relatively high sensitivity for radial distances 
as low as 4-6 fm, well inside the strong absorption radius. 
Deeper inside this radial range, the refractive index, de- 
fined as n = 



= Jl- 



is almost real and reaches values 



as high as n = 2.6, comparable to that of diamond. 

As mentioned already, it was shown in Refs. H, 0, 0| 
that the elastic scattering of light heavy ion systems such 
as ^^0-|-^^C and ^^0-1- ^^O shows sufficient transparency 
for the cross section to be dominated by far-side scat- 
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FIG. 2: (Color online) Notch test for the radial sensitivity of 
the potential. The bare Xo is calculated for the first potential 
in Tableim A Gaussian spike with a 10% strength and a width 
(FWHM) of 0.2 fm is superimposed on the bare potential at 
various radii. Points are connected by lines to guide the eye. 



tering. Structures appearing in the elastic scattering 
angular distributions at intermediate angles have been 
identified as Airy minima of a nuclear rainbow, due to 
a destructive interference between two far-side trajecto- 
ries which sample the interior of the potential. At 19 
MeV/nucleon the ^Li scattering data show rapid, diffrac- 
tive Fraunhofer oscillations at small angles due to the 
strong near- far amplitude interference (Fig. |2J|. Beyond 
the crossover the near-side amplitude makes a negligi- 
ble contribution to the cross section. The shoulder and 
the deep minimum seen at 9 McV/nucleon (Fig. ^ are 
washed out in the far-side amplitude and only a broad, 
less pronounced minimum survives, followed by a broad 
Airy maximum and an exponential, structureless decay 
of the cross section at large angles. Clearly, both the data 
at 9 and 19 MeV/nucleon (Figs. ^andl^J show far-side 
dominance as a possible manifestation of refractive ef- 
fects. However, this simple dominance does not explain, 
by itself, the difference in the angular distributions seen 
at these energies, suggesting a difference in the reaction 
mechanism. In fact the above picture has been already 
challenged by Anni 12] and by Michel et al. p3| for the 
simple reason that the far-side amplitude has never been 
decomposed in subamplitudes which would explain the 
quoted interference. We come back to this topic in Sec- 
tion ^T] For the moment we adopt the interpretation of 
Michel et al. [s^l and denote the complex structure at in- 
termediate angles in our data as pre-rainbow oscillations. 

In the remainder of this section we discuss the ability 
of the folding model to describe the pre-rainbow oscilla- 
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FIG. 3: (Color online) Woods-Saxon optical model analysis 
of elastic data at 130 MeV ( Table UTI ) . Far-side/near-side 
cross sections are also shown by dashed and dotted lines, re- 
spectively. The lower curves are multiplied by 10~^ and 10"'', 
respectively. 

tion seen at 9 MeV/nucleon and the rainbow patterns at 
higher energies. Data at somewhat lower energies are also 
examined in order to see if the plateau feature persists in 
adjacent systems and on a larger energy range. Our pre- 
ferred model is the nuclear matter approach of Jeukenne, 
Lejeune and Mahaux (JLM) which incorporates a 
complex, energy and density dependent parametrization 
of the NN effective interaction obtained in a Brueck- 
ner Hartree-Fock approximation from the Reid soft core 
nucleon-nucleon potential. The systematic study Q of 
the elastic scattering between p-shell nuclei at energies 
around 10 MeV/nucleon led to the surprising result that 
on average, the imaginary part of the folded JLM poten- 
tial was perfectly adequate to describe such reactions and 
did not need any renormalization {Nw = 1.00 ± 0.09), 
while the real component needed a strong renormaliza- 
tion, in line with other effective interactions used in fold- 
ing models. However, the present data extend to a much 
larger angular range and need further refinements of this 
model. 

In the JLM model the complex form factor for the 
optical potential is given by 

C/(i?)= j dr^dT2pM)p2{r2)v{p,E,s)g{s) (3) 

where v is the (complex) NN interaction, Pi(2) ^-re the 
single particle densities of the interacting partners, cal- 
culated in a standard spherical Hartree-Fock procedure 
using the energy density functional of Beiner and Lom- 
bard with the surface term adjusted to reproduce the 
total binding energy [sil [s^ , s = ri + R — r2 \s the NN 



separation distance between interacting nucleons and p 
is the overlap density. The effective NN interaction con- 
tains an isovcctor component which gives a negligibly 
small contribution for p-shell nuclei but is included here 
for convenience in conjunction with appropriate single 
particle isovector densities. The smearing function g{s) 
is taken as a normalized Gaussian 0, Il9l l33| , 

f(^) = ^3^exp(-sVi') (4) 

which tends to a (5-function for t — s- 0, while for finite 
values of the range parameter t it increases the rms ra- 
dius of the folding form factor by = (3/2)t^, leaving 
unchanged the volume integral. Inclusion of a smear- 
ing function with a varying range parameter, greatly in- 
creases the ability of the folding form factor to simulate 
the radial dependence of DPP. 

The geometric or arithmetic mean of the overlapping 
densities has been used to define the overlap density p in 
Eq. El 

P=[Pi{ri + \^P2{r2-\^Y'^ (5) 

and 

P = + ^s) +P2(^2 - ^s)]- (6) 

The former was introduced by Campi and Sprung in 
density-dependent Hartree-Fock calculations (3^. It is 
physically appealing since the overlap density tends to 
zero when one of the interacting nucleons is far from the 
bulk, and to the nuclear matter saturation value at com- 
plete overlap. The approximation in Eq. (Sjis similar to 
that used in folding calculations with density-dependent 
M3Y effective interactions (35j|, except for the factor 1/2 
which has been introduced here because JLM interac- 
tion is defined only up to the nuclear matter saturation 
value p < Pq- It has been suggested to us 36j that 
the drawbacks seen in our earlier analysis of the scat- 
tering of ''Li at 19 MeV/nucleon (see Fig. 6a in Rcf. 
3]) may be due to the weak density dependence intro- 
duced by Eq. El and thus rainbow patterns could not 
be reproduced. However, the optical model analysis pre- 
sented above showed clearly that the pre-rainbow oscil- 
lations (at 9 MeV/nucleon) and rainbow patterns (at 19 
MeV/nucleon) could be described if and only if the po- 
tentials have the proper rms radius. It turns out that the 
smearing procedure described above is essential in simu- 
lating the complicated radial dependence of the dynamic 
polarization potential. 

In the earlier analysis Q], fixed values for the range 
parameters ty — 1.2 fm and tyy — 1.75 fm, found from 
a global analysis of the data were used. Only the renor- 
malization factors Nv and Nw were left free in the fits 
for each case. In the present analysis with double folded 
potentials, all four parameters: two strength parameters 
{Ny and Nw) and two range parameters {ty and tw), 
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FIG. 4: (Color online) Comparison of the JLMl folding model calculations (full lines) with present data at 9 MeV/nucleon. 
The parameters are given in Table Hill Far-side (dashed) and near-side (dotted) cross sections are indicated in the form of 
ratios to the Rutherford cross sections. 



have been searched simultaneously to fit the data for each 
case 

Udf (r) = NvV(r, ty) + iN^Wiv, t^) (7) 

to obtain a phenomenological representation of the DPP 
as a uniform renormalization of the depths and radii of 
the folding potentials. The calculations using Eqs. © 
and l|njl are dubbed JLMl and JLM2 respectively. As 
these give very similar results only JLMl parameters are 
listed in Table UTTl and the results of the calculations are 
shown in Figs. |3|to[7| At 9 MeV/nucleon (Fig. the 
same pattern emerges as with Woods- Saxon form fac- 
tors. The pre-rainbow oscillation is carried entirely by 
the dominant far-side component. Some other high or- 
der structures appear at angles near 180° as the result of 
near/far amplitude interference. At most forward angles 



this interference produces an inner Fraunhofer crossing 
which give rise to a deep minimum in the cross section. 

For ^Li-t-^Be at 63 MeV, JLMl calculation failed to 
describe the oscillation near 9 — 80° for the simple rea- 
son that data required a rms radius for the real potential 
Rv = 3.4 fm, while the bare JLM interaction predicts a 
minimal i?v=3.6 fm for tn « 0. This once again reflects 
the critical role played by the radial behavior of DPP. 
This is also illustrated in the upper right quadrant of Fig. 
|S1 where two JLM solutions for the reaction -'^^B-l-^Be at 
10 MeV/nucleon are indicated (see also Table ^Oj. The 
solution with smaller real volume integral which better 
fits the forward angles predicts a smooth, exponentially 
decaying cross section beyond 9 « 60° . The second solu- 
tion with a real volume integral close to the critical value 
Jvcrit ~ 300 MeVfm^^ gives rise to a shallow pre-rainbow 
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FIG. 5: (Color online) Comparison of JLMl folding model calculations with ^Li scattering data at 19 MeV/nucleon. The data 
are taken from 3]. Two JLMl solutions are indicated for ^^B+^Be reaction. The parameters are given in Table Far-side 
(dashed) and near-side (dotted) cross sections are indicated in ratio to Rutherford cross sections. 



oscillation at these angles (not covered by experiments). 
The high selectivity of the pre-rainbow oscillations to the 
optical potentials is also illustrated in Fig. Elwhere other 
^Li scattering data from literature, at somewhat lower 
energies, are explored. The ^Li-t-^^C data at 50 MeV 
(isf could be described in the whole angular range only 
with potentials exceeding the critical value of the real vol- 
ume integral found before. In Fig. [7|we show ^Li-|-^^C 
elastic scattering data at 7 energies between 15 and 50 
MeV/nucleon. Now, even at high energy (Fig. [7^ the 
JLMl description of the rainbow patterns is exemplary 
(to be compared with Figs. 6 a) and c) of Ref. Q). This 
suggests that the geometrical details of the optical po- 
tential rather than the density dependence are essential 
for a correct description of ^■''Li elastic scattering at low 
and intermediate energies. 



A close examination of the parameters in Table UTTl re- 
veals an erratic variation of the range parameters iy(vi/) 
from one energy to another and from system to system. 
As mentioned above, this largely reflects the mass and en- 
ergy dependence of DPP. The other parameters are more 
stable. The strength parameter Ny decreases slowly from 
5 to 16.5 MeV/nucleon and then increases again up to 53 
MeV/nucleon, the highest energy at which reliable data 
exist. This may suggest that DPP reaches its maximum 
amplitude at energies around 16 MeV/nucleon. On aver- 
age the Ny values in Table ITTll are somewhat larger than 
in our earlier analysis |3| reflecting the need for stronger 
refractive effects, but again Nyy approaches unity, on av- 
erage. 
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FIG. 6: (Color online) Comparison of the JLMl folding model calculations with ®Li scattering data on light targets at 30 and 
50 MeV laboratory energy. For the ^^C target at 50 MeV, only a solution with a real volume integral exceeding the critical 
value Jv=SOO MeV fm"^ fTable llTTt is able to reproduce both forward and intermediate angles (right bottom panel). Far-side 
(dashed lines) and near-side (dotted) cross sections are indicated in ratio to Rutherford cross sections. 



IV. TRANSFER REACTION 

As already mentioned, in one experimental run 
we have also measured the neutron transfer reaction 
"C(^Li,8Li)i2c at E(^Li)=63 MeV. The purpose of the 
study was to determine the ANC for the ground state of 
^Li, and then, using charge symmetry to relate it with 
that in its mirror nucleus ^B. The ANC was then used to 
calculate the astrophysical factor 6*17 that gives the rate 
of the proton capture reaction ^Be(p, 7)^B, of crucial im- 
portance for the solar neutrino problem. The major ad- 
vantage of the neutron transfer reaction over its mirror 
proton transfer reaction is that it involves a stable beam, 
and, therefore, a much more precise and detailed angular 
distribution could be measured. That allowed the deter- 



mination of the admixture of the minor component lpi/2 
in the wave function of the ground state of *Li (and ^B, 
respectively), dominated by the IP3/2 orbital. The re- 
sults of this experiment were reported in Ref. I^^. In 
that study we paid particular attention to the depen- 
dence of the results on the optical model potentials used 
in the entrance and exit channels. 

Eleven different combinations of entrance/exit poten- 
tials were used to show that the resulting values for C^^^^ 
and Cp_^^^ are very stable, when the potentials are rea- 
sonable. The potentials used were either volume Woods- 
Saxon forms with the parameters from similar projectile- 
target combinations at similar energies, or were obtained 
from the double folding procedure with the renormaliza- 
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FIG. 7: (Color online) Comparison of JLMl folding model 
calculations (full lines) with high energy ^Li scattering data. 
The source of data and calculation parameters are given in 
Table |TTT1 The far-side (dashed lines) cross sections are in- 
dicated in ratio to the Rutherford cross sections (and from 
the top curve below, each case is multiplied by an extra 10^^ 
factor). The near-side cross section (not shown) is important 
only at most forward angles. 



tion coefficients from the previous paper Calcula- 
tions done after the publication with the new (deeper) 
potential " 227" of Table IIIII in both entrance and exit 
channels lead to minor (~5%) variations in the results. 
The very good agreement between the experimental data 
and the DWBA calculations and between the results of 
present and previous calculations (Figure |SJ) shows that 
the region of the potential contributing to transfer (the 
surface) is well and uniquely described. This simulta- 
neous description of elastic and transfer data is also an 
argument for the complete determination of the optical 
potentials. 



V. DISPERSION RELATION 

The dispersion relation is a fundamental property of 
the optical potential (see for example [42]) and a se- 
lection between ambiguous potentials can be performed 
by studying the dispersive properties of these poten- 
tials, provided accurate analyses of experimental data are 
available over a large energy range. 

The threshold anomaly which manifests itself as a 
sharp increase of the real optical potential for energies 
close to the Coulomb barrier, has been explained by Na- 
garajan, Mahaux and Satchler as due to the opening 
of reaction channels with increasing energy. An applica- 
tion of the dispersion relation for elastic scattering of ^^O 
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FIG. 8: (Color online) The angular distributions for the neu- 
tron transfer reaction ^■^C(^Li,*Li)^^C to the ground (top, 
open points) and first excited state (bottom, full points) of 
*Li. The calculations shown (full line) are done using the po- 
tential "227" in Table El The data are shown as points and 
the separate contributions of the pi/2 — > P3/2 (dashed line)and 
Pi/2 P3/2 (dotted line) components are shown in both cases. 



on ^"^Pb at energies around 80 MeV accounted well for 
this effect. Later it was conjectured by Mahaux, Ngo and 
Satchler that for loosely bound nuclei, this anomaly 
may be absent. For these nuclei, the strong coupling with 
breakup channels gives rise to a repulsive DPP which 
compensates the strong attractive component. Accord- 
ing to Sakuragi ^] this effect would explain the large 
renormalization needed by most of the effective inter- 
actions used in the folding model for elastic scattering 
of ^'^Li. The coupling with inelastic channels alone has 
been invoked by Gomez-Camacho et al. fi^ to explain 
this reduction. The earlier analysis of Kailas (i^l found 
strong dispersive effects for ^Li-|-^^C scattering. Recent 
studies by Tiede et al. and by Pakou et al. of 
^Li+^^Si at near barrier energies found that the strength 
of the real part of the folding potential using the M3Y 
interaction remains almost independent of energy, sug- 
gesting a cancellation between the attractive (dispersive) 
component and the strong repulsive dynamic polarization 
potential arising from the coupling to continuum states. 
Another study of ^^'^U+'^^^Vh near the Coulomb barrier 
^3 found that at low energies the DPP is of opposite sign 
for the two projectiles and there is a threshold anomaly 
for ""Li but none for ^Li. No significant fusion hindrance 
caused by breakup effects was found in the fusion reac- 
tion of *^'^Li on a Co target near the Coulomb barrier 
[48l |. thus leading the authors to conclude that breakup 
suppression above the barrier appears to be a common 
feature of ^'^Li induced reactions. 
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Therefore, the energy dependence of the ^^^Li optical 
potential is far from clear and the competition between 
dispersive (attractive) and coupling to continuum (re- 
pulsive) effects need to be studied more carefully. An 
earlier study showed that the total reaction cross 
section for ^Li scattering saturates at energies around 
20 MeV/nucleon and therefore dispersive effects could 
be identified by accumulating good optical potentials in 
this energy range. The real and imaginary volume in- 
tegrals for the optical potentials obtained in the previ- 
ous sections are plotted in Fig. O Both Woods-Saxon 
and JLM folding results have been included. These are 
supplemented with values derived from the smooth OMl 
potential of Trcka et al. • 

We assume that the local optical potential may be writ- 
ten as = Vq+AV{E) where Vq is independent of energy 
and AV{E) is the energy dependent DPP. We ignore the 
spurious energy dependence of Vq arising from non local- 
ity which is expected to be weak for heavy ions. We use 
the dispersion relation connecting the imaginary and real 
volume integrals in the subtracted form, 



Jav,eSE) ^ {E ~ Es)'^ I — 



Jw{E') 



{E' - Es){E' - E) 



dE' (8) 



where Es is a reference energy and V is the principal 
value of the integral. In principle the evaluation of this 
equation requires the knowledge of Jw values at all en- 
ergies. The above subtracted form takes advantage of 
the fact that the energy dependence of Jw far from sat- 
uration energy is not very important and the unknown 
contributions are absorbed by normalizing to the empir- 
ical value at a convenient reference energy. 



Jav^eAE) = Jav{E) - Jav{Es) 



(9) 



Two schematic models have been employed here to es- 
timate the energy dependence of the imaginary volume 
integral. A first one approximates this energy depen- 
dence by straight line segments which makes the 
evaluation of Eq. [S] analytical. A more realistic energy 
dependence is given by 



JwiE) ^ J^il - l3eM-aE)) 



(10) 



where the parameters Jy^^^lTO MeV fm^, a=0.023 
MeV~^ and (3—0.95 describe better the energy depen- 
dence in the important range 0-20 MeV/nucleon. In both 
calculations the reference energy was set to £^s=156 MeV, 
an energy where the JLM folding model gives precise 
values for volume integrals. In general, the calculated 
dispersion contributions get more repulsive as the en- 
ergy increases, and the corresponding real potentials get 
shallower, in qualitative consistency with phenomenol- 
ogy. An empirical logarithmic dependence of the form 
Jy = -785 + 951n(i;) has been found in Ref. |50| mostly 
based on unique OM potentials determined from 35 and 
53 MeV/nucleon ^Li scattering on light targets. This 
matches perfectly the dependence obtained with the dis- 
persion relation for E>10 MeV/nucleon, but disagrees 



at lower energies. In fact, this logarithmic dependence 
is physically meaningful and can be understood on the 
basis of the dispersion relation with a schematic (line 
segments) approach for the imaginary volume integral. 
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FIG. 9: (Color online) Energy dependence of the real (solid 
points) and imaginary (open points) volume integrals ob- 
tained in the analyses with Woods-Saxon and folding (JLM) 
optical potentials. The stars show the values obtained from 
the OMl optical potential of Ref. jl^. The curves for Jy 
are the result from the dispersion relation, normalized to the 
empirical value at 26 AMeV, assuming the schematic mod- 
els shown for Jw- The dash-dotted line gives the empirical 
energy dependence of the real volume integral of Ref. 
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A relatively strong localized energy variation is 
predicted by the linear model in the range 0-20 
MeV/nucleon, while the exponential model predicts a 
smooth dependence on the entire range of energies. This 
last calculation is much closer to the data and seem to 
confirm Jy=320 MeV fm'^ as the most realistic value at 9 
MeV/nucleon, in surprising agreement with values found 
for the more bound system ^^0-1-^^0 (see e.g. Fig 6 in 
reL |H3)- Most probably the phenomenological values 
found at 5 MeV/nucleon are due to the erratic variation 
in the WS parameters due to the rapidly changing elastic 
scattering angular distributions 52] near the resonance 
energy region around 20 MeV. 



VI. SEMICLASSICAL BARRIER AND 
INTERNAL BARRIER AMPLITUDES 



Once we have established the main features of the av- 
erage OM potential, we turn now to study the reaction 
mechanism in the elastic scattering of ^'^Li on light tar- 
gets at 9 MeV/nucleon using semiclassical methods. The 
far-side dominance observed in the angular distributions 
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at 9 and 19 MeV/nucleon is not able to explain the differ- 
ences in the reaction mechanism at these energies. The 
reason is of course that the far/near (F/N) decomposition 
method does not perform a dynamic decomposition of 
the scattering function, but merely decomposes the scat- 
tering amplitude into traveling waves. The intermediate 
angle structures, such as those observed in our angular 
distributions, have been repeatedly interpreted as arising 
from the interference of two ranges in angular momenta 
^< and £> contributing to the same negative deflection 
angle. However, the corresponding cross sections (Jf< 
and aF> cannot be isolated because their dynamic con- 
tent (S'-matrix) is not accessible. 

The semiclassical uniform approximation for the scat- 
tering amplitude of Brink and Takigawa is well 
adapted to describe situations in which the scattering 
is controlled by at most three active, isolated, complex 
turning points. An approximate multireflection series ex- 
pansion of the scattering function can be obtained, the 
terms of which have the same simple physical meaning as 
in the exact Debye expansion for the scattering of light 
on a spherical well. The major interest in this theory 
comes from the fact that it can give precious informa- 
tion on the response of a nuclear system to the nuclear 
interior. Recent application of this technique helped 
to clarify the controversial problem of " Airy oscillation" 
seen in low energy scattering (sj. 
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where the Langer prescription has been used for the cen- 
trifugal term. This guarantees the correct behavior of 
the semiclassical wave function at the origin Then 
we calculate the deflection function. 



^W = --^ 2 I ,r (12) 

where r\ is the outer zero of the square root, i.e. the 
radius of closest approach to the scatterer and \x is the 
reduced mass. Note that with the replacement = 
by/2iiE, Eq. 1121 becomes identical with the classical de- 
flection function Q{b), where b is the impact parame- 
ter. The results are shown in Fig. 1101 The behavior 
of 6(A) is the one expected for a strong nuclear poten- 
tial in a near orbiting kinematical situation in which the 
cm. energy approximately equals the top of the bar- 
rier for some specific angular momentum. The deflec- 
tion functions exhibit no genuine minimum, but rather 
a pronounced cusp close to an orbiting logarithmic sin- 
gularity. Therefore any interpretation of structures in 
angular distributions in terms of Airy oscillations can be 
discarded. Rather we need an interpretation appropri- 
ate for orbiting, a well documented situation in classical 
physics (54| . We identify the cusp angular momenta as 
orbiting momenta (Aq) since they are related with the 
coalescence of two (barrier) turning points and the in- 
nermost turning point given by the centrifugal barrier 
become classically accessible. There are two branches 
that can be distinguished, an internal branch, for low 
active momenta A < Ao related to semiclassical trajecto- 
ries which penetrate into the nuclear pocket and a less 
developed external (barrier) branch (A > Ao) related to 
trajectories deflected at the diffuse edge of the potential. 
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FIG. 10: (Color online) Semiclassical deflection functions for 
the real potentials shown in Table |n] In the top left quad, 
glories are indicated by open dots and the orbiting angular 
momentum by an arrow. 

For the potentials in Tabled (the flrst entry for each of 
the four cases measured here) we discard the absorptive 
terms and define the effective potential as, 

K//(r)=nO + ^^, A = £+i (11) 



However this simple calculation cannot determine the 
relative importance of these branches and provides no 
information about the interference effects of the cor- 
responding semiclassical trajectories. To clarify these 
points it is best to go into the complex r-plane and look 
for complex turning points, i.e. the complex roots of the 
quantity Ec.m. — Veff — iW. This is an intricate numer- 
ical problem, because, for a WS optical potential, the 
turning points are located near the potential singulari- 
ties and there are an infinite number of such poles. The 
situation for integer angular momenta is depicted in Fig. 
^for the reaction ^Li-|--^^C at 54 MeV using the poten- 
tial "225" in Table HH Only turning points nearest the 
real axis are retained and we observe an ideal situation 
with three, well isolated turning points for each partial 
wave. Even small absorption plays an essential role in 
the motion of turning points. Removing the imaginary 
part W, the barrier turning points (ri,2) become complex 
conjugates while the internal turning point is purely real 
(open symbols in Fig. [TT|l . 

The multireflection expansion of the scattering func- 
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FIG. 11: (Color online) Complex turning points (full symbols) 
for the potential "225" shown in Table HTl at integer angular 
momenta. Open symbols denote turning points for the real 
potential alone. Stars indicate complex poles of the potential. 



tion in the Brink- Takigawa approach reads, 



q=0 



(13) 



where, 



exp(2i(5^) 



and for q ^ 0, 



N{S2i/n) 
1 exp [ 2i{qS32 + S21 + 6i) 



(14) 



(15) 



In these equations is the WKB (complex) phase shift 
corresponding to the turning point ri, N{z) is the barrier 
penetrabihty factor. 



Niz) 



r{z- 



J— exp {zhiz — z) 



(16) 



and Sij is the action integral calculated between turning 
points ri and r^, 

S^J - £ dr{^[E,.,^. - K// - iW\Yl^ (17) 

and 5*32 are independent of the integration path pro- 
vided they lie on the first Riemann sheet and collision 
with potential poles is avoided. Each term in Eq. 1131 
has a simple physical interpretation. The first term (the 



barrier term, denoted also Sb) retains contributions from 
trajectories reflected at the barrier, not penetrating the 
internal region. The gth term corresponds to trajectories 
refracted q times in the nuclear interior with q-\ reflec- 
tions at the barrier turning point r2 ■ Summation of terms 
g > 1 can be recast into a single term. 



exv\2i{S32 + ^21 + ^i)] 



1 



A^(52i/7r) 



1 



•exp[2z532]/^(52iA) 
(18) 

and is known as the internal barrier scattering function. 
When the absorption in the nuclear interior is large, 
the second factor in the above equation reduces to one 
and we are left with the expression used in |30l |. Since 
the semiclassical scattering function is decomposed addi- 
tively, SwKB = Sb + Si, the corresponding total scatter- 
ing amplitude is decomposed likewise as fwKB = Jb + Ji 
and conveniently the corresponding barrier and internal 
barrier angular distributions are calculated as asj — 
using the usual angular momentum expansion of 
the amplitudes. 
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FIG. 12: (Color online) Semiclassical decomposition of sca- 
tering function for the WS potential of Fig. 1111 Barrier (open 
circles) and internal barrier components (squares) are indi- 
cated. The exact total quantum S-matrix is indicated by 
small dots. The line is a cubic spline interpolation of the 
total semiclassical scattering function for the same potential. 

The accuracy of the semiclassical calculation has been 
checked by comparing the barrier and internal barrier ab- 
sorption profiles with the exact quantum-mechanical re- 
sult in Fig. El First, one observe that the semiclassical 
B /I expansion is an exact decomposition of the quantum 
result. They are virtually identical at the scale of the 
figure. The internal component gets significant values up 
to the grazing angular momentum (€g = 15) and is neg- 
ligibly small beyond this value. The barrier component 
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resembles a strong absorption profile and this justifies the 
interpretation that it corresponds to that part of the flux 
not penetrating into the nuclear interior. For values near 
the orbiting angular momentum (£o ~12), the two com- 
ponents interfere and a downward spike appears in the 
total profile, in complete agreement with the quantum 
result. Second, the B/I components are almost decou- 
pled in the angular momentum space and therefore they 
will contribute in different angular ranges. 

Semiclassical cross sections are compared with the data 
in Figinifor the reaction ^Li+^'^C at 54 MeV. Better in- 
sight into this technique is obtained by further decom- 
posing the B/I components into far and near (BF/BN 
and IF/IN) subcomponents. Clearly, the barrier com- 
ponent dominates the forward angle region. Fraunhofer 
diffractive oscillations appear as the result of BF and BN 
interference. At large angles, the internal contribution 
accounts for the full cross section. As both B/I contribu- 
tions are dominated by the far-side component (Fig. 1131 
bottom panels), we show in Fig. E|the angles at which 
the phase difference of the BF and IF amplitudes passes 
through an odd multiple of tt, i.e. where minima should 
be expected. Since the crossing angle (where ctb « cr/) is 
about « 75° and lies just in between predicted minima, 
the coherent interference around this angle gives rise to 
the " plateau" (constructive) and the deep minimum (de- 
structive) at 9 Ki 80° . Similar consideration apply to the 
other three reactions. 

Thus, the intermediate angle exotic structure in an- 
gular distributions for the elastic scattering of ^'^Li on 
light targets can be understood as a result of coherent 
interference of two far-side subamplitudes generated by 
different terms in the uniform multireflection expansion 
of the scattering function (terms q=0 and q—l in Eq. 
I13|l . corresponding to the scattering at the barrier and 
the internal barrier. This interference effect appears as 
a signature of a surprisingly transparent interaction po- 
tential for loosely bound nuclei ®'^Li which allows part 
of the incident flux to penetrate the nuclear interior and 
reemerge with significant probability. 



VII. CONCLUSIONS 

We have performed precise measurements on extended 
angular ranges of the elastic scattering of loosely bound 
nuclei ^'^Li on I2,i3(^ g^j-^j 9gg f^^j. projectile-target 
combinations at 9 MeV/nucleon and reanalyzed previous 
data for the scattering of ^Li at 19 MeV/nucleon in an ef- 
fort to obtain systematic information on the interaction 
of p-shell nuclei with light targets. Optical potentials 
for these nuclei are needed for studies in which highly 
peripheral transfer reactions involving radioactive nuclei 
are used as indirect methods for nuclear astrophysics and 
are an important factor in the accuracy and reliability of 
these methods. At the present status of the experimental 
techniques, the best information on the optical potentials 
for radioactive nuclei can be obtained only by extrapo- 



lation from adjacent less exotic nuclei. Our intention is 
to narrow the ambiguities in the optical model potentials 
by systematic studies of the scattering of loosely bound 
projectiles over a large range of angles and energies, and 
extract information that can be used for systems involv- 
ing radioactive projectiles, for which elastic scattering 
data of very good quality are not easily available. We 
demonstrate this procedure by reanalyzing the one neu- 
tron transfer reaction ^^C(^Li,^Li)^^C using optical po- 
tentials obtained in the present study. 

The present data, which extend over a much larger an- 
gular range than previously measured, confirm the exis- 
tence of an exotic intermediate angle structure, observed 
pre viously by Trcka et al. It was interpreted in Ref. 
|18l | as a diffractive effect arising from an angular mo- 
mentum dependent absorption. We adopt an opposite 
point of view and interpret these structures as refrac- 
tive effects arising from a fine balance between the real 
and imaginary components of the optical potential. We 
have performed a traditional analysis of our data in terms 
of Woods-Saxon and microscopic JLM folded potentials. 
Both approaches lead to the conclusion that the opti- 
cal potential is deep and surprisingly transparent, in line 
with findings for more bound systems. Folding model 
form factors have been renormalized in the usual way in 
order to account for the energy and radial dependence of 
the dynamic polarization potential. It is suggested that 
DPP attains its maximum amplitude at approximately 
16 MeV/nucleon for these systems. The intermediate 
angle structures could be reproduced only with poten- 
tials exceeding a critical volume integral of about 300 
MeV fm'^ and, consequently, are severely selective, lim- 
iting the ambiguities in the determination of the OMP. 
The remaining discrete ambiguities could be removed by 
a dispersion relation analysis. Based on a good estima- 
tion of the absorption at low energy (5-20 MeV/nucleon), 
this analysis allowed us to extract a smooth energy de- 
pendence of the optical potential. Our analysis did not 
find any spectacular anomaly near the Coulomb barrier 
and seems to confirm, to some extent, the conjecture of a 
canceling effect between the repulsive dynamic polariza- 
tion potential due to the coupling with breakup channels 
and the attractive, dispersive component of the optical 
potential. 

In our previous study Q we found a simple recipe 
to obtain OMP for loosely bound p-shell nuclei from a 
double folding procedure using the JLM effective NN 
interaction. The already independent real and imag- 
inary parts were smeared with constant, but different 
ranges ty = 1.2 fm and tyy = 1.75 fm, which ac- 
counted for the well known need for a wider imaginary 
potential to describe the experimental data. We found 
that a considerable renormalization of the real part was 
needed Ny = 0.37 ± 0.02 (leading to volume integrals 
Jy ~ 220 MeV fm'^), but not for the imaginary part 
Nw = 1.00 ± 0.09. That recipe was already successfully 
applied to predict the elastic scattering angular distribu- 
tions of RNBs on light targets in a number of cases at 
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FIG. 13: (Color online) Semiclassical barrier and internal barrier decomposition of the cross section. The turning points and 
scattering function are those from Figs. Illl and lT^ respectively. Each component is further dcomposed into far-side (dashed) 
and near-side (dotted) components. 



energies around 10 MeV/nucleon. The present analysis 
shows that in order to reproduce the structures observed 
at intermediate angles in the same cases measured, one 
needs to allow for a more complicated radial dependence 
of the dynamic polarization potential, energy and target 
dependent, and require deep real potentials with volume 
integrals larger than a critical value of Jvcrit ~ 300 MeV 
fm'' . This is a conclusion of the phenomenological analy- 
ses and is supported by the dispersion relation analysis. 
However, the elastic scattering data in the angular range 
of the Fraunhofer oscillations and the transfer reactions 
can be equally well described by the previous potentials 
produced by the folding procedure with fixed smearing 
ranges for the effective NN interaction and the simple 
renormalization of Ref. , showing that the potentials 
are well described in the surface region. 



In an effort to clarify the reaction mechanism respon- 
sible for the intermediate angle structures found at 9 
MeV/nucleon, we performed extensive semiclassical cal- 
culations within the uniform multireflection expansion of 
the scattering function of Brink and Takigawa. It has 
been shown that using complex trajectories, the (exter- 
nal) barrier/internal barrier expansion is an exact realiza- 
tion of the dynamic decomposition of the quantum result 
into components responsible for that part of the incident 
flux reflected at the barrier and the part of the flux which 
penetrates into the nuclear interior and reemerges with 
significant probability. By combining the B/I decom- 
position with the usual far-side/near-side expansion, we 
explain the intermediate angle structure as a coherent 
interference effect of two subamplitudes (BF and IF). 
Thus, this refractive effect appears as a signature of a 
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FIG. 14: (Color online) The phase difference of the far-side 
barrier/far-side internal barrier amplitudes as a function of 
scattering angle. Large dots indicate the predicted interfer- 
ence minima. For easier comparison, the experimental cross 
section is shown as 7 -|- \og{a/aR) to match the scale of the 
figure. 



highly transparent interaction potential. 
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TABLE I: List of the elastic scattering experiments presented 
in this paper. 



No. 


Reaction 


E [MeV] 0iab[deg.] 


1 


'^Li + '■'C 


54 


2 - 56 


2 


"Li + 


54 


2 - 59 


3 


^Li + ^Be 


63 


4 - 52 


4 


^Li + 


63 


4 - 56 


5 


^Li + ^Be 


130 


4-47 


6 


^Li + 


130 


4 - 47 



TABLE II: Best fit Woods-Saxon parameters. Reduced radii arc defined in the heavy ion convention. All lengths arc given in 
fin, depths and energies in MeV, cross sections in mb and volume integrals in MeV fm^. Coulomb reduced radius is fixed to 
rc=l fm. Rv and Rw are the rms radii of the real and imaginary potentials, respectively. 



Reaction 


Energy 


M) 


Wo 


rv 


rw 0.V aw 


x' 




Jv 


Rv 


Jw 


Rw 




[MeV] 


[MeV] 


[MeV] 


[fm] 


[fm] [fm] [fm] 




[mb] [MeV fm^] 


[fm] 


[MeV fm^] 


[fm] 


"Li+i'C 


54 


225.47 


15.75 


0.503 


1.157 0.900 0.737 


17.71 


1309 


338 


3.70 


121 


4.59 






371.31 


17.70 


0.439 


1.109 0.856 0.777 


13.60 


1322 


419 


3.47 


125 


4.56 


^Li+i^C 


54 


225.28 


14.75 


0.502 


1.181 0.916 0.707 14.62 


1327 


327 


3.76 


114 


4.63 






364.46 


16.95 


0.443 


1.133 0.871 0.744 


14.24 


1338 


403 


3.53 


119 


4.58 


^Li+«Be 


63 


225.85 


24.74 


0.536 


0.941 0.828 0.980 


10.14 


1456 


369 


3.49 


146 


4.66 






368.34 


29.38 


0.478 


0.882 0.790 1.004 


11.85 


1470 


464 


3.28 


153 


4.62 


^Li+i^C 


63 


227.94 


15.37 


0.529 


1.186 0.932 0.669 


20.09 


1367 


328 


3.87 


107 


4.64 






278.86 


24.19 


0.594 


1.050 0.789 0.721 


20.04 


1334 


411 


3.53 


126 


4.38 


^Li+i3c 


130 


149.11 


29.73 


0.636 


0.932 0.885 0.929 


2.61" 


1403 


282 


3.90 


132 


4.62 


■^Li+^Be 


130 


143.41 


33.64 


0.581 


0.829 0.892 1.094 


3.03" 


1446 


295 


3.76 


169 


4.80 


i^C+^Be 


130 


159.85 


24.43 


0.674 


0.983 0.868 0.914 


13.69 


1552 


280 


3.96 


104 


4.79 



" uniform 10% errors. 



19 



TABLE III: Best fit JLMl parameters. The notations are those from the text. Lengths are given in fm, energies in MeV, cross 
sections in mb and volume integrals in MeV fm'^. 



Reaction 


Energy 


tv 


tw 


AT 


7\r 

Nw 


X 




Jv 




7 

Jw 


r> 

Uw 




[Me V \ 


[fmj 


[fmJ 








[mb] [MeV fm^'] 


[fmJ 


[MeV im ] 


[fmJ 


LiH- 


30" 


0.30 


2.45 


0.60 


0.46 


14.8 


1371 


396 


3.66 


72 


4.93 




50'' 


0.08 


2.78 


0.56 


0.78 


12.0 


1315 


373 


3.64 


120 


4.42 




54 


0.08 


2.76 


0.54 


0.77 


21.4 


1556 


351 


3.64 


116 


5.11 




90*" 


0.70 


2.70 


0.52 


1.24 


18.4 


1591 


313 


3.73 


173 


4.96 




99^^ 


0.60 


1.75 


0.47 


1.01 


4.21 


1225 


277 


3.69 


145 


4.27 




124^* 


0.60 


1.75 


0.51 


1.09 


3.96 


1243 


292 


3.69 


168 


4.28 




15o 


0.50 


1.50 


0.50 


0.94 


7.98 


1146 


271 


3.66 


154 


4.19 






0.60 


1.75 


0.58 


1.11 


5.87 


1231 


305 


3.68 


185 


4.28 




210-^ 


0.20 


1.35 


0.56 


0.93 


23.5 


1062 


276 


3.59 


161 


4.05 




318'' 


0.80 


1.95 


0.60 


0.85 


9.00 


1069 


251 


3.69 


148 


4.35 


^Li+"C 


54 


0.08 


2.76 


0.54 


0.77 


21.4 


1556 


351 


3.64 


116 


5.11 


^Li+^'^O 


50'' 


0.50 


2.81 


0.55 


0.60 


13.4 


1643 


346 


3.64 


91 


5.24 


^Li+^Be 


63 


0.09 


1.20 


0.46 


0.98 


19.5 


1538 


274 


3.64 


152 


4.80 


^Li+«C 


63 


0.12 


2.59 


0.52 


0.78 


19.0 


1652 


335 


3.74 


113 


5.07 


^Li+i3c 


130' 


0.13 


1.97 


0.48 


1.02 


4.58 


1392 


280 


3.73 


146 


4.50 


^Li+^Be 


130' 


0.12 


2.34 


0.50 


1.23 


7.98 


1404 


304 


3.62 


183 


4.65 


14N+13C 


162 


1.44 


1.82 


0.39 


0.73 


33.1 


1563 


220 


4.29 


89 


4.66 


i"B+9Be 


100 


1.89 


1.02 


0.30 


1.01 


6.9 


1266 


185 


4.33 


146 


4.08 






0.47 


2.28 


0.48 


0.93 


29.6 


1558 


298 


3.75 


133 


4.79 



data from ' (sl, " g^, = [Z^/ [H, ^ [H, '' [H, ' uniform 10 % errors. 



